arXiv:1508.03700v2 [cond-mat.stat-mech] 9 Mar 2016 


Extreme event statistics of daily rainfall: Dynamical systems approach 


G. Cigdem Yalcin 

Department of Physics, Istanbul University, 34134, Vezneciler, Istanbul, Turkey 
Pau Rabassa and Christian Beck 

Queen Mary University of London, School of Mathematical Sciences, Mile End Road, London El 4^8, UK 

We analyse the probability densities of daily rainfall amounts at a variety of locations on the 
Earth. The observed distributions of the amount of rainfall fit well to a (/-exponential distribution 
with exponent q close to (/ « 1.3. We discuss possible reasons for the emergence of this power 
law. On the contrary, the waiting time distribution between rainy days is observed to follow a 
near-exponential distribution. A careful investigation shows that a (/-exponential with q « 1.05 
yields actually the best fit of the data. A Poisson process where the rate fluctuates slightly in a 
superstatistical way is discussed as a possible model for this. We discuss the extreme value statistics 
for extreme daily rainfall, which can potentially lead to flooding. This is described by Frechet 
distributions as the corresponding distributions of the amount of daily rainfall decay with a power 
law. On the other hand, looking at extreme event statistics of waiting times between rainy days 
(leading to droughts for very long dry periods) we obtain from the observed near-exponential decay of 
waiting times an extreme event statistics close to Gumbel distributions. We discuss superstatistical 
dynamical systems as simple models in this context. 


I. INTRODUCTION 

The statistical analysis of precipitation data is an in¬ 
teresting problem of major environmental importance [H- 
13 . Of particular interest are extreme events of rainfall, 
which can lead to floodings if a given threshold is ex¬ 
ceeded. From a mathematical and statistical point of 
view, it is natural to apply extreme value statistics to 
measured rainfall data, but it is not so clear which class 
of the known extreme value statistics, if any, is applica¬ 
ble, and how the results differ from one spatial location 
to another. Another interesting quantity to look at is 
the waiting time between rainy days. Extreme events of 
waiting times in this context correspond to droughts. So 
an interesting question is what type of drought statistics 
is implied by the observed distribution of waiting times 
between rainfall periods if this is extrapolated to very 
long waiting times. 

In this paper we present a systematic investigation of 
the probability distributions of the amount of daily rain¬ 
fall at variety of different locations on the Earth, and of 
waiting time distributions on a scale of days and hours. 
Our main result is that the observed distributions of 
the amount of rainfall are well-fitted by (/-exponentials, 
rather than exponentials, which suggests that techniques 
borrowed from nonextensive statistical mechanics @ and 
superstatistics 0 could be potentially useful to better 
understand the rainfall statistics. An entropic exponent 
of (/ Ri 1.3 is typically observed. In fact, based on the fact 
that (/-exponentials asymptotically decay with a power 
law, we discuss the corresponding extreme value statis¬ 
tics that is highly relevant in the context of floodings 
produced by extreme rainfall events. We also investigate 
the waiting time distribution between rainy events, which 
is much better described by an exponential, although an 
entropic exponent close to 1 such as q ~ 1.05 seems to 
give the best fits. We discuss possible dynamical reasons 


for the occurrence of (/-exponentials in this context. 

One possible reason could be superstatistical fluctua¬ 
tions of a variance parameter or a rate parameter. Let 
us explain this a bit further. Superstatistical techniques 
have been discussed in many papers [Ml and they rep¬ 
resent a powerful method to model and/or analyse com¬ 
plex systems with two (or more) clearly separated time 
scales in the dynamics. The basic idea is to consider for 
the theoretical modelling a superposition of many sys¬ 
tems of statistical mechanics in local equilibrium, each 
with its own local variance parameter /3, and finally per¬ 
form an average over the fluctuating (3. The probability 
density of /3 is denoted by /(/3). Most generally, the pa¬ 
rameter P can be any system parameter that exhibits 
large-scale fluctuations, such as energy dissipation in a 
turbulent flow, or volatility in financial markets. Another 
possibility is to regard P as the rate parameter of a local 
Poisson process, as done, for example, in (2^. Ultimately 
all expectation values relevant for the complex system un¬ 
der consideration are averaged over the distribution f{P). 
Many applications have been described in the past, in- 
cluding modelling the statistics of classical turbulent flow 
quantum turbulen ce 1^ . space-time granular¬ 
ity 12^ ■ stock price changes [l^jwind velocity fluctua¬ 
tions [^, sea level fluctuations [^, infection pathways 
of a virus [ 2 ^ , migration trajectories of tumour cells , 
and much more pdl. . Superstatistical systems, 

when integrated over the fluctuating parameter, are ef¬ 
fectively described by more general versions of statistical 
mechanics, where formally the Boltzmann-Gibbs ent ropy 
is replaced by more general entropy measures [l^, Il7j |. 
The concept can also be generalized to general dynami¬ 
cal systems with slowly varying system parameters, see 
for some recent rigorous results in this direction. 

Our main goal in this paper is to better understand the 
extreme event statistics of rainfall at various example lo¬ 
cations on the Earth. We will start with a careful analysis 
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of experimentally recorded time series of the amount of 
rainfall measured at a given location, whose probability 
density is highly relevant to model the corresponding ex¬ 
treme event statistics [ssl - fs^ . Ultimately of course all 
this rainfall dynamics can be formally regarded as be¬ 
ing produced by a highly nonlinear and high-dimensional 
deterministic dynamical system in a chaotic state, pro¬ 
ducing the occasional rainfall event, hence it is useful to 
keep in mind the basic results of extreme event statis¬ 
tics for weakly correlated events as generated by mix¬ 
ing dynamical systems. Recently there has been much 
activity on the rigorous application of extreme values 
theory to deterministic dynamical systems [iol - li^ and 
also to stochastically perturbed ones [i§-[5lj. A remark¬ 
able feature of the dynamical system approach is that 
there exist some correlations between events, and hence 
the extreme value theory used to tackle it must account 
for this correlation going beyond a theory that is just 
based on sequences of events that are statistically inde¬ 
pendent. In the superstatistics approach, correlations are 
also present, due to the fact that parameter changes take 
place on long time scales, but the relaxation time of the 
system is short as compared to the time scale of these 
parameter changes, so that local equilibrium is quickly 
reached. 

What is worked out in this paper is a comparison with 
experimentally measured rainfall data, to decide which 
extreme event statistics should be most plausibly applied 
to various questions related to amount of rainfall and 
waiting times between rainfall events. Extreme value the¬ 
ory quite generally tells us (under suitable asymptotic in¬ 
dependence assumptions) that there are only three possi¬ 
ble limit distributions, namely, the Gumbel, Frechet and 
Weibull distribution. But are these assumptions of near¬ 
independence satisfied for rainfall data, and if yes, which 
of the above three classes are relevant? This is the subject 
of this paper. We will also discuss simple deterministic 
dynamical system models that generate superstatistical 
processes in this context. 

The paper is organized as follows. In Section II we 
present histograms of rainfall statistics, extracted from 
experimentally measured time series of rainfall at various 
locations on the Earth. What is seen is that the prob¬ 
ability density of amount of rainfall is very well fitted 
by q-exponentials. We discuss the generalized statistical 
mechanics foundations of this based on nonextensive sta¬ 
tistical mechanics with entropic index q, with q « 1.3. 
In section III we look at waiting time distributions (on 
a daily and hourly scale) between rainy episodes. These 
are observed to be close to exponential functions, similar 
as for the Poisson process. However, a careful analy¬ 
sis shows that a slightly better fit is again given by q- 
exponentials, but this time with q much closer to 1. A 
simple superstatistical model for this is discussed in sec¬ 
tion IV, a Poisson process that has a rate parameter that 
fluctuates in a superstatistical way. We review standard 
extreme event statistics in section V and then, in sec¬ 
tion VI, based on the measured experimental results of 


rainfall statistics, we develop the corresponding extreme 
value statistics. In section VII we analyse the ambigui¬ 
ties that arise for the extreme event statistics of waiting 
times, depending on whether we assume the waiting time 
distribution is either an exact exponential or a slightly 
deformed ^-exponential as produced by superstatistical 
fluctuations. Finally, in section VIII we describe a dy¬ 
namical systems approach to superstatistics. 

II. DAILY RAINFALL AMOUNT 

DISTRIBUTIONS AT VARIOUS LOCATIONS 

We performed a systematic investigation of time se¬ 
ries of rainfall data for 8 different example locations on 
the Earth (Figs. 1-8). The data are from various pub¬ 
licly available web sites. When doing a histogram of the 
amount of daily rainfall observed, a surprising feature 
arises. All distributions are power law rather than expo¬ 
nential. They are well fitted by so-called g-exponentials, 
functions of the form 

eg{x) = {1 + {1 - q)x)'^^^^~‘>'> ( 1 ) 

which are well-motivated by generalized versions of sta¬ 
tistical mechanics relevant for systems with long-range 
interactions, temperature fluctuations and multifractal 
phase space structure Of course the ordinary ex¬ 

ponential is recovered for g —5> 1. Whereas the data of 
most locations are well fitted by g « 1.3, Central England 
and Vancouver have somewhat lower values of g closer to 
1.13. 

One may speculate what the reason for this power law 
is. Nevertheless, the formalism of nonextensive statis¬ 
tical mechanics @ is designed to describe complex sys¬ 
tems with spatial or temporal long-range interactions, 
and g-exponentials occur in this formalism as generalized 
canonical distributions that maximize g-entropy 

= ( 2 ) 

^ i 

where the pi are the probabilities of the microstates i. 
Ordinary statistical mechanics is recovered in the limit 
g —>■ 1, where the g-entropy Sq reduces to the Shannon 
entropy 

Si = -'^Pi\ogPi. (3) 

i 

The generalized canonical distributions maximize the g- 
entropy subject to suitable constraints. In our case the 
constraint is given by the average amount of daily rain¬ 
fall at a given location. The way rainfall is produced 
is indeed influenced by highly complex weather systems 
and condensation processes in clouds, so one may spec¬ 
ulate that more general versions of statistical mechanics 
could be relevant as an effective description. Also for hy¬ 
drodynamic turbulent systems [3, IMl and pattern form¬ 
ing systems [111 these generalized statistical mechanics 



Central England - (1931 - 2014) 



FIG. 1: Observed probability distribution of amount of daily 
rainfall (minus average value) in Central England (red data 
points). The blue curve is a fit with a g-exponential eq{—x). 
The parameter q is fitted for the various locations. For Cen¬ 
tral England, q = 1.13. 


Darwin - (1975- 2014) 



FIG. 2: As Fig. 1, but for Darwin {q — 1.30). 


methods have previously been shown to yield a good ef¬ 
fective description. The amount of rain falling on a given 
day is a complicated spatio-temporal stochastic process 
with intrinsic correlations, as rainy weather often has a 
tendency to persist for a while, both spatially and tem¬ 
porally. The actual value of q for the observed rainfall 
statistics reflects characteristic effective properties in the 
climate and temporal precipitation pattern at the given 
location. For temperature distributions at the same lo¬ 
cations as in Fig.1-8, see (^ . 


III. WAITING TIME DISTRIBUTIONS 
BETWEEN RAINY DAYS 

Another interesting observable that we extracted from 
the data is the waiting time distribution between rainy 
episodes. We did this both for a time scale of days and 


Dijbai - (1968 - 2014) 



FIG. 3: As Fig. 1, but for Dubai {q — 1.30). 


Eureka - (1951 - 2014) 



FIG. 4: As Fig. 1, but for Eureka {q = 1.30). 


a time scale of hours. A given day is marked as rainy 
if it rains for some time during that day. The waiting 
time is then the number of days one has to wait until 
it rains again, this is a random variable with a given 
distribution which we can extract from the data. Results 
for the waiting time distributions are shown in Fig. 9-14. 


Hongkong - (1997 - 2014) 



FIG. 5: As Fig. 1, but for Hongkong {q — 1.30). 
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Ottawa - (1939 - 2014) 



FIG. 6: As Fig. 1, but for Ottawa (q = 1.15). 


Sydney - (1910 - 2014) 



FIG. 7: As Fig. 1, but for Sydney (q = 1.25). 


Central England - waiting times 



FIG. 9: Waiting time distribution between rainy days for Gen- 
tral England (red data points). The blue curve is a fit with a 
q-exponential. For waiting time distributions, the parameter 

q is Tnnrb rlnspr to 1 fn = 1 O.'i for Dpntral Envlfinrl) 


Darwin - waiting times 



What one observes here is that the distribution is nearly 
exponential. That means the Poisson process of nearly 
independent point events of rainy days is a reasonably 
good model. 

At closer inspection, however, one sees that again a 
slightly deformed q-exponential, this time with q Ri 1.05, 
is a better fit of the waiting time distribution. As worked 


FIG. 10: As Fig. 9, but for Darwin (q = 1.05). 

out in the next section, one may explain this with a super- 
statistical Poisson process, i.e. a Poisson process whose 
rate parameter - on a long time scale- exhibits fluctua¬ 
tions that are y^-distributed, with a rather large number 
of degrees of freedom n. 


Vancouver - (1937 - 2014) 



FIG. 8: As Fig. 1, but for Vancouver (q = 1.13). 


IV. SUPERSTATISTICAL POISSON PROCESS 

We start with a very simple model for the return time 
of rainfall events (or extreme rainfall events) on any given 
time scale. This is to assume that the events follows a 
Poisson process. For a Poisson process the waiting times 
are exponentially distributed, 

p(t|/3) = /3exp(-/3t). 

Here, t is the time from one event (peak over threshold) 
to the next one, and ,5 is a positive parameter, the rate 
of the Poisson process. The symbol p{t\f3) denotes the 
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Eureka - waiting times Vancouver - waiting times 




FIG. 11: As Fie. 9. but for Eureka Co = LlOF 


FIG. 14: As Fig. 9, but for Vancouver {q = 1.10). 


Ottawa - waiting times 



FIG. 12: As Fig. 9, but for Ottawa (q = 1.05). 


conditional probability density to observe a return time 
t provided the parameter /3 has certain given value. 

The key idea of the superstatistics approach SHI can 
be applied to this simple model, thus constructing a su- 
perstatistical Poisson process. In this case the parameter 
/3 is regarded as a fluctuating random variable as well, but 
these fluctuations take place on a very large time scale. 
For example, for our rainfall statistics the time scale on 


Sydney — waiting times 



FIG. 13: As Fig. 9, but for Sydney (q — 1.06). 


which /3 fluctuates may correspond to weeks (different 
weather conditions) whereas our data base records rain¬ 
fall events on an hourly basis. 

If /3 is distributed with probability density /(/3), and 
fluctuates on a large time scale, then one obtains the 
marginal distribution of the return time statistics as 


POO POO 

p{t)= f{l3)p{t\l3)= /(/3)/3exp(-^t). (4) 

Jo Jo 


This marginal distribution is actually what is recorded 
when we sample histograms of the observational data. 

By inferring directly on a simple model for the distribu¬ 
tion /(/3), a more complex model for the return times can 
be derived without much technical complexity. For exam¬ 
ple, consider that there are n different Gaussian random 
variables Xi,i = 1 ,..., n, that influence the dynamics of 
the intensity parameter /3 as a random variable. We may 
thus assume as a very simple model that /3 = 
with E{Xi) = 0 and E{Xf) ^ 0. Then the probability 
density of P is given by a y^-distribution: 


fm 


1 

r(n/2) 



nil 

exp 



where n is the number of degrees of freedom and is a 
shape parameter that has the physical meaning of being 
the average of /3 formed with the distribution /(/3). 

The integral (|3]) is easily evaluated and one obtains the 
g-exponential distribution: 

p(t)^(l + 6(g-l)t)i/(i-«), 


where g = 1 -1- 2/(n -I- 2) and b = — g). 

To sum up, this model generates g-exponential distri¬ 
butions by a simple mechanism, fluctuations of a rate 
parameter fj. Typical g-values obtained in our fits are 
g = 1.3 for rainfall amount and g = 1.05 for waiting time 
between rainfall events. 
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V. EXTREME VALUE THEORY FOR 
STATIONARY PROCESSES 

Classic extreme value theory is concerned with the 
probability distribution of unlikely events. Given a sta¬ 
tionary stochastic process Xi,X 2 ,..., consider the ran¬ 
dom variable defined as the maximum over the first 
n-observations: 

M„ = max(Xi,... ,X„). (5) 

In many cases the limit of the random variable may 
degenerate when n —^ oo. Analogously to central limit 
laws for partial sums, the degeneracy of the limit can be 
avoided by considering a rescaled sequence a„(M„ — 5„) 
for suitable normalising values a„ > 0 and G M. In¬ 
deed, extreme value theory studies the existence of nor¬ 
malising values such that 

P{an{M„-bn) <x) ^ G{x). (6) 

as n —>■ oo, with G{x) a non-degenerate probability dis¬ 
tribution. 

Two cornerstones in Extreme Value Theory are the 
Fisher-Tippet Theorem [sj and the Gnedenko Theorem 
[55l| . The former asserts that if the limiting distribu¬ 
tion G exist, then it must be either one of three pos¬ 
sible types, whereas the latter theorem gives necessary 
and sufficient conditions for the convergence of each of 
the types. A third cornerstone in Extreme Value The¬ 
ory are the Leadbetter conditions [H, . These are a 

kind of weak asymptotic independence conditions, under 
which the two previous theorems generalize to station¬ 
ary stochastic series satisfying them. Let us review these 
results in somewhat more detail. 

In the case where the process Xi is independent iden¬ 
tically distributed (i.i.d.) the Fisher-Tippett Theorem 
states that if Ai, A2 ,... is i.i.d. and there exist sequences 
> 0 and € M such that the limit distribution G is 
non-degenerate, then it belongs to one of the following 
types: 

Type I : G{x) = exp (—e““) for cc G K.. This distribu¬ 
tion is known as the Gumbel extreme value distri¬ 
bution (e.v.d.). 

Type II : G{x) = exp(—for a: > 0; G{x) = 0, 
otherwise; where a > 0 is a parameter. This family 
of distributions is known as the Frechet e.v.d. 

Type III: G{x) = exp (—(—a;)“), for x < 0; G{x) = 1, 
otherwise; where a > 0 is a parameter. This family 
is known as the Weibull e.v.d. 

A further extension of this result is the Gnedenko The¬ 
orem, which provides a characterization of the conver¬ 
gence in each of these cases. Let Ai, A2 ,... be an i.i.d. 
stochastic process and let F be its cumulative distribu¬ 
tion function. Gonsider xm = sup{x| F'(x) < 1}. The 


following conditions are necessary and sufficient for the 
convergence to each type of e.v.d.: 

Type I : There exists some strictly positive function h{t) 
such that lim, - = e~^ for all real x; 

Type II: xm = +00 and limt^oo = X with 

a > 0, for each x > 0; 

Type III: xm < 00 and limt^o = a:“, with 

a > 0, for each x > 0. 

This result implies that the extremal type is completely 
determined by the tail behaviour of the distribution F{x). 

VI. EXTREME EVENT STATISTICS FOR 
EXPONENTIAL AND Q-EXPONENTIAL 
DISTRIBUTIONS 

The rainfall data were well described by q- 
exponentials, but waiting time distributions were ob¬ 
served to be close to ordinary exponentials, with q only 
deviating by a small amount from 1. Let us now discuss 
the differences in extreme value statistics that arise from 
theses different distributions. 

In the case where A is distributed as an ordinary ex¬ 
ponential function with parameter A, we have 

exp(—At) if t > 0 . . 

1 if t<0. 

It is not difficult to check that the exponential distri¬ 
bution belongs to the Gumbel domain of attraction. In 
other words, the extreme events associated to the expo¬ 
nential distribution will be Gumbel distributed. 

Recall that the g-exponential function is defined as 

exp^(f) := [1 -b (I - , 

with 1 < g < 2. A random variable A is g-exponential 
distributed (with parameter A) if its density function is 
equal to (2 — g)Aexp^(—Ax). In such a case, its hazard 
function is 

1 _ / exp,. i-Xt/q') if t > 0 , . 

1 if t<0. 

where g' = 1/(2 — g). 

Using the Gnedenko theorem if follows that the g- 
exponential distribution belongs to the Frechet domain 
of attraction. In this case the shape parameter of the 
Frechet distribution a is equal to ^Ef- 

VII. MODEL UNCERTAINTY 

Extremely large waiting times for rainfall events corre¬ 
spond to droughts. Clearly, it is interesting to extrapo¬ 
late our observed waiting time distributions to very large 


1 - F{t) = I 
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FIG. 15: Horizontal axis: Multiple k of the mean, vertical 
axis: probability P for an event at that level. 

time scales. However, in Section IlIIl we saw that in most 
cases it is difficult to discern if the waiting time distri¬ 
bution is that of a Poisson process, distributed as an 
exponential, or if it is a ^-exponential with q close to 
one. This can make a huge difference for extreme value 
statistics. The aim of this section is to assess the impact 
of choosing one or the other model. 

Consider k a constant and X a random variable mod¬ 
elled either by an exponential or a ^-exponential. To nor¬ 
malize the problem, we can scale our analysis in terms of 
the mean. In other words, we look at the probability of 
X being bigger than a multiple, say fc-times, the mean of 
X. 

If X is distributed like an exponential with parameter 
A, its mean is equal to 1/A and its hazard function is 
given by ©. Then it is easy to check that 

P{X > kE{X)) = exp (—fc). 

On the other hand, if X is distributed like a q- 
exponential with parameter A, its mean is equal to 
1/A(3 — 2q) (provided 1 < q < 3/2) and its hazard func¬ 
tion is given by (|S]). In this case we have 

P{X > kE{X)) = expg, . 

Recall that the exponential distribution can be under¬ 
stood as the limit of the q-exponential as q goes to 1. 
This is also true for the probability above, which con¬ 
verges to exp(—fc) as q goes to 1. In Fig 15 we plotted 
the probability of an event of level k for different values 
of q. For instance the probability of having an observa¬ 
tion bigger than 5 times the mean is 0.0068 for q = 1, 
0.0095 when q = 1.05 and 0.026 when q = 1.3. When we 
look at the probability of an observation bigger than 10 
the mean, it is 0.000045, 0.00023 and 0.0068 respectively. 
Apparently, the predicted drought statistics is very dif¬ 
ferent choosing either the value 9 = 1 or the very similar 
value q = 1.05 for the observed waiting time distribution. 


This illustrates the general uncertainty in model building 
for extreme rainfall and drought events [5^ . 

VIII. DYNAMICAL SYSTEMS APPROACH 

Ultimately, the weather and rainfall events at a given 
location can be regarded as being produced by a very 
high-dimensional deterministic dynamical system ex¬ 
hibiting chaotic properties. It is therefore useful to ex¬ 
tend the superstatistics concept to general dynamical sys¬ 
tems, following similar lines of arguments as in the recent 
paper Q. 

The basic idea here is that one has a given dynamics 
(which, for simplicity, we take to be a discrete mapping 
/a : U —>• U on some phase space U) which depends on 
some control parameter a. If a is changing on a large time 
scale, much larger than the local relaxation time deter¬ 
mined by the Perron-Frobenius operator of the mapping, 
then this dynamical system with slowly changing control 
parameter will ultimately generate a superposition of in¬ 
variant densities for the given parameter a. Similarly, 
if we can calculate return times to certain particular re¬ 
gions of the phase space for a given parameter a, then in 
the long term the return time distribution will have to 
be formed by taking an average over the slowly varying 
parameter a. Clearly, the connection to the previous sec¬ 
tions is that a rainfall event corresponds to the trajectory 
of the dynamical system being in a particular subregion 
of the phase space f2, and the control parameter a corre¬ 
sponds to the parameter /3 used in the previous sections. 

Let us consider families of maps fa depending on a con¬ 
trol parameter a. These can be a priori arbitrary maps 
in arbitrary dimensions, but it is useful to restrict the 
analysis to mixing maps and assume that an absolutely 
continuous invariant density Pa(x) exists for each value 
of the control parameter a. The local dynamics is 

Xn+l = faiXn)- (9) 

We allow for a time dependence of a and study the long¬ 
term behavior of iterates given by 

Xn = fa„0 fa„_^0 ...fa^{xo)■ (10) 

Clearly, the problem now requires the specification of the 
sequence of control parameters ai,..., a„ as well, at least 
in a statistical sense. One possibility is a periodic orbit 
of control parameters of length L. Another possibility 
is to regard the Oj as random variables and to specify 
the properties of the corresponding stochastic process in 
parameter space. This then leads to a distribution of 
parameters a. 

In general, rapidly fluctuating parameters Oj will lead 
to a very complicated dynamics. However, there is a sig¬ 
nificant simplification if the parameters Oj change slowly. 
This is the analogue of the slowly varying temperature 
parameters in the superstatistical treatment of nonequi- 
librium statistical mechanics [^. The basic assumption 




of superstatistics is that an environmental control pa¬ 
rameter a changes only very slowly, much slower than 
the local relaxation time of the dynamics. For maps this 
means that significant changes of a occur only over a 
large number T of iterations. For practical purposes one 
can model this superstatistical case as follows: One keeps 
oi constant for T iterations (T >> 1), then switches af¬ 
ter T iterations to a new value 02, after T iterations one 
switches to the next values 03, and so on. 

One of the simplest examples is a period-2 orbit in the 
parameter space. That is, we have an alternating se¬ 
quence 01,02 that repeats itself, with switching between 
the two possible values taking place after T iterations. 
This case was given particular attention in [s^, and rig¬ 
orous results were derived for special types of maps where 
the invariant density Pa as a function of the parameter o 
is under full control, so-called Blaschke products. 

Here we discuss two important examples, which are of 
importance in the context of the current paper, namely 
how to generate (in a suitable limit) a superstatistical 
Langevin process, as well as a superstatistical Poisson 
process, using strongly mixing maps. 

Example 1 Superstatistical Langevin-like pro¬ 
cess We take for fa a map of linear Langevin type 
[H, [5^. This means fa is a 2-dimensional map given 
by a skew product of the form 

Xn+l = g{Xn) (11) 

yn+i = - 5 ) ( 12 ) 

Here g denote the average of iterates of g. It has been 
shown in that for r —5> 0, t = nr finite this deter¬ 
ministic chaotic map generates a dynamics equivalent to 
a linear Langevin equation [^ . provided the map g has 
the so-called (p-mixing property and regarding the 
initial values Xq G [0,1] as a smoothly distributed ran¬ 
dom variable. Consequently, in this limit the variable yn 
converges to the Ornstein-Uhlenbeck process [H, and 
its stationary density is given by 



The variance parameter /3 of this Gaussian depends on 
the map g and the damping constant a. If the parameter 
a changes on a very large time scale, much larger than the 
local relaxation time to equilibrium, one expects for the 
long-term distribution of iterates a mixture of Gaussian 
distributions with different variances /3~^. For example, 
a period 2 orbit of parameter changes yields a mixture of 
two Gaussians 



Generally, for more complicated parameter changes on 
the long time scale T, the long-term distribution of iter¬ 
ates pn will be given by a mixture of Gaussians with a 


suitable weight function /i(/3) for r —> 0: 

p(y) ^ J df3 (15) 

This is just the usual form of superstatistics used in 
statistical mechanics, based on a mixture of Gaussians 
with fluctuating variance with a given weight function 
0. Thus for this example of skew products the super¬ 
statistics of the map fa reproduces the concept of super¬ 
statistics in nonequilibrium statistical mechanics, based 
on the Langevin equation. In fact, the map fa can be 
regarded as a possible microscopic dynamics underlying 
the Langevin equation. The random forces pushing the 
particle left and right are in this case generated by de¬ 
terministic chaotic map g governing the dynamics of the 
variable Xn- Generally it is possible to consider any ip- 
mixing map here [s^ . Based on functional limit theo¬ 
rems, one can prove equivalence with the Langevin equa¬ 
tion in the limit r ^ 0. 

We note in passing that if the mixing property is not 
satisfied, then of course much more complicated proba¬ 
bility distributions are expected for sums of iterates of 
a given map g. Numerical evidence has been presented 
that in this case often g-Gaussians with g > 1 are ob¬ 
served 

Example 2 Superstatistical Poisson-like process 

Take /„ : H ^ H to be a strongly mixing map, for 
example the binary shift map fa{x) = 2x mod 1 on the 
interval H = [0,1]. Gonsider a very small subset of the 
phase space of size e, say 7^ = [1 — e, 1] and a generic 
trajectory of the binary shift map for a generic initial 
value Xq = o’i2“b where ai G {0,1} are the dig¬ 

its of the binary expansion of the initial value Xq. We 
define a ‘rainfall’ event to happen if this trajectory en¬ 
ters /e (of course, for true rainfall events the dynamical 
system is much more complicated and lives on a much 
higher-dimensional phase space). It is obvious that the 
above sequence of events follows Poisson-like statistics, 
as the iterates of the binary shift map are strongly mix¬ 
ing, which means asymptotic statistical independence for 
a large number of iterations. Indeed between successive 
visits of the very small interval 1 ^, there is a large number 
of iterations and hence near-independence. Hence the bi¬ 
nary shift map generates a very good approximation of 
the Poisson process for small enough e, and the waiting 
time distribution between events is exponential. 

We may of course also look at a more complicated sys¬ 
tem fa, where we iterate a strongly mixing map fa which 
depends on a parameter a, and where the invariant den¬ 
sity Pa of the map fa depends on the control parameter in 
a nontrivial way. Examples are Blaschke products, stud¬ 
ied in detail in [^ . If the parameter a varies on a large 
time scale, so does the probability Pe = /j Pa{x)dx of 
iterates to enter the region 7^, and hence the rate param¬ 
eter of the above Poisson-like process will also vary. The 
result is a superstatistical Poisson-like process, generated 
by a family of deterministic chaotic mappings fa- In this 
way we can build up a formal mathematical framework to 
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dynamically generate superstatistical Poisson processes. 

IX. CONCLUSION 

We started this paper with experimental observations: 
The probability densities of daily rainfall amounts at a 
variety of locations on the Earth are not Gaussian or ex¬ 
ponentially distributed, but follow an asymptotic power 
law, the g-exponential distribution. The corresponding 
entropic exponent q is close to q ~ 1.3. The waiting 
time distribution between rainy episodes is observed to 
be close to an exponential distribution, but again a care¬ 
ful analysis shows that a g-exponential is a better fit, 
this time with g close to 1.05. We discussed the corre¬ 
sponding extreme value distributions, leading to Gumbel 
and Frechet distributions. We made contact with a very 
important concept that is borrowed from nonequilibrium 
statistical mechanics, the superstatistics approach, and 


pointed out how to generalize this concept to strongly 
mixing mappings that can generate Langevin-like and 
Poisson-like processes for which the corresponding vari¬ 
ance or rate parameter fluctuates in a superstatistical 
way. Of course rainfall is ultimately described by a very 
high dimensional and complicated spatially extended dy¬ 
namical system, but simple model systems as discussed 
in this paper may help. 
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